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ABSTRACT 

Wiggins, Edwin George, M.S . in Nuclear Engineering, Purdue 
University , January 1969. An Investigation of the Use of the Hot-Wire 
Anemometer in Non-Isothermal Air Flow . Major Professor: Alexander 

Ses onske . 

Procedures were developed for the calibration of the hot-wire 
anemometer in a non- isothermal flow. Various methods of using the re- 
sults of calibration to compute turbulent fluctuations of velocity and 
temperature were considered, and the Arya and Plate modification of the 
Kovasznay Fluctuation Diagram method was found to give the most accurate 
results. The method requires independent measurement of the temperature 
fluctuations which was obtained with a rapid response thermocouple. 

The calculated values of velocity turbulence intensity show devi- 
ations of 36 percent from the results of Laufer; however, this deviation 
is to be expected in view of the estimated experimental error. These 
errors were probably caused by contamination of the wire during the 
lengthy calibration process. 

All measurements were made in air flowing in a 1.0625 inch diameter 
pipe at Reynolds Numbers between 30,000 and 50,000. A heat input of 430 
watts over the last 58 inches before the probe was used to generate 
temperature differences between wall and centerline of 48 to 63^ at the 
probe location. 

The procedures developed and the experimental problems identified 
are the contribution of this investigation to a comprehensive research 
program of non- isotherma 1 hot-wire measurements in fluids of a wide range 



of Prandtl Numbers. 
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INTRODUCTION 

The hot-wire anemometer has been used extensively for the measure- 
ment of velocity fluctuations in isothermal flow of gases. The principles 

I 2 3 

and procedures involved are discussed by Kovasznay ’ and by Kronauer . 

Corrsin has shown that in addition, the hot-wire anemometer has the 
potential to measure simultaneous velocity and temperature fluctuations 

in non- isothermal flow. Experimental work of this nature has been done 

5 6 7 8 

by Corrsin and Uberoi , Gibson Chen and Lin , Deissler , Kunstman , and 

9 

Arya and Plate . The present work was concerned with the development of 
experimental procedures and the identification of problems associated with 
making of such measurements in a non- isotherma 1 air flow. 

In isothermal flow the only flow variable affecting the anemometer 
output is the velocity. Thus the anemometer voltage fluctuation is 
directly proportional to the velocity fluctuation and the anemometer 
steady voltage is directly proportional to the mean velocity. 

In non- isothermal flow, however, both the velocity and ambient 
temperature fluctuations affect the anemometer output. Considerable dif- 
ficulty has been encountered by previous workers in separating the two 

g 

effects. Errors up to 750 percent have been reported in non- isotherma 1 
measurements, while in isothermal measurements 10 percent error or less 
is expected. 

The non- isothermal flow measurements of the present work were pre- 
ceeded by velocity fluctuation measurements in isothermal flow with a 



2 



hot-wire paralleling the work of Laufer and Sandborn and 
fluctuation measurements in non- isotherma 1 flow with a rapid 
thermocouple. Thus the velocity and temperature fluctuation 
measured with the hot-wire in non- is otherma 1 flow could be c 
results of more direct measurements. 



temperature 
response 
intens it ies 
ompared with 



THEORETICAL BACKGROUND 



The constant temperature hot-wire anemometer consists of a fine 
tungsten or platinum wire maintained by a feedback system at some constant 
temperature above that of the fluid in which it is inserted. The con- 
vective heat loss from the wire is balanced by joulean heating of the 

wire by the current supplied by the feedback system. 

12 

King solved the convection problem by assuming potential flow 

around the wire and obtained the following relationship; 

1 T -T ^ 

Nu = (A + BRe"^) (I + h - ) , 

a 

where Re = Pup . 

Although it is recognized that King's assumptions are unacceptable, the 
functional form of his relationship has proved to be a satisfactory re- 
presentation of experimental data. Thus it is seen that the joulean 
heating current (or voltage) required to maintain the wire at a constant 
temperature is a function of the instantaneous velocity of the flow and 
the difference in temperature between wire and fluid. 

In isothermal flow, the temperature difference is constant and the 
system can be used to measure the mean and fluctuating components of the 
instantaneous velocity. In non- is othermal flow the fluid temperature 
has a fluctuating component that affects the heat transfer in addition 



to the effect of the fluctuating component of the velocity. Thus in 



4 



principle, the system can be used to measure velocity and temperature 
fluctuations, since it is sensitive to both fluctuations. 



The 

hot-wire 

Nu 



Hot-Wire in Isothermal Flow 

following parallels Corrsin*s^ treatment for a constant 
anemometer. For heat transfer to flow over a cylinder: 

(A + B [1 + — - ] , (1) 

a 



current 



In the range of interest, 



1 

2 




1 . 



( 2 ) 



Thus Equation 1 may be simplified to 



Nu 



A + B Re 



1/2 



(3) 



By definition , 

(A) 



( 5 ) 

Substitution of Equations 4 and 5 into 3 yields 

2 

R (t" - T ) ■ (6) 

w w a ^ ^ 

By definition , 

Re = ^ . (7) 

P 



Nu = 



hD 





k(T - T ) 
w a 



For electrical resistance heating, 

o"- 

^ R ttD£ ■ 

w 



5 



Thus , 



oD 1/2 

+ (TTkiiB [ ^ ] ) ^ 



w w 



( 8 ) 



The quantities in parentheses on the right side of Equation 8 are 
constants. For constant temperature operation of the hot-wire in an 
isothermal flow the following mean and deviation quantities are used: 



u = u -f u ’ 



and 



E = E + e’ 



Substitution of Equations 9 into Equation 8 yields, after small 
terms have been discarded » 



72 . , 1/2 

) = (nkilA) + (-nk^B [ ^ 1 ) u (1 + -^ ) 

w^ w a 



By definition, 



(9) 



( 10 ) 



T - T = — (R - R ) . 
w a aR w a 
o 



( 11 ) 



According to the binomial theorem 

. 1/2 , , 

(1 + — ) = 1 + ^ ^ + . . . . 

u 2 u 



( 12 ) 



Since << 1, it is assumed that higher order terms may be 

neglected in Equation 12. Thus Equation 10 becomes 









2Ee’ 



,TTk£ 



1/2 



R (R -R ) R (R -R ) "aR 
w w a w w a o 



= (^ A) + (21^ B[^] )(u^''^+-^). 



aR ‘ M 
o 



.-1/2 

2u 



(13) 



Time averaging of Equation 13 yields 






R (R -R ) 
w w a 



. ( IM A) + ( iM E 



aR 



aR 



1/2 

pD -1/2 

— ] )u . 



(1^) 
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Subtraction of Equation 14 from 13 gives 



2£e^ 

R (R - R ) 
w w a 



ctR p 



1/2 



- 1/2 

2u 



or 

4Ee^ 

— 2^/2 — * 

R (R ~ R ) Cu u 

w w a 

Note that 

CD 

uR \i 

o 

— 

and that this same quantity appears as the coefficient of u^ in 
Equation 14. Squaring, time averaging, and taking the square 
root of Equation 16 yields 

4¥/ e’^ ^ u' ^ 

— 1/2 — 

R (R - R ) Cu ' u 

w w a 

Rapid Response Thermocouple in Non- Isothermal Flow 
The behavior of a thermocouple over a fairly large range of 
temperature may be expressed as 

E = E + aT . 
o 

The following fluctuating quantities are introduced: 

E = E + e’ and T = T + t* . 

Substitution of Equations 20 into 19 yields 

E+e’=E +a(T+t’). 
o 



(15) 



(16) 



(17) 



(18) 



(19) 



( 20 ) 



( 21 ) 
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Time averaging of Equation 21 yields 

E = E + a? . 
o 



( 22 ) 



Subtracting Equation 22 from 21 yields 
e’ = at* 



(23) 



Squaring, time averaging, and taking the square root of Equation 23 
gives ■ mZ 

/ e’2 = a/ t'2 . (24) 



Hot-Wire in Non-1 sothermal Flow 

4 

The following parallels Corrsin’s treatment for a constant 
current hot wire anemometer. 

Equation 8 is still valid, however, P and k are no longer constants. 
Substitution of Equation 11 into 8 yields 



R (R -R ) 
w w a 



= (• 



7T£ 

aR 



A)k +(■ 



tt£D 



1/2 



aR 



-C M 
P 






(2f ) 



It is assumed that 

k ^ 

C y Pr 
P 

is constant, and that C itself is also constant. 

P 

The following mean and deviation quantities are considered: 



Ra = + Ta’ P = P + p’ 

k=k+k E=E+e' 

u = u + u ' 



( 26 ) 
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Substitution of Equations 26 into 2S yields 

2 



(E+e') 



R (R -R -r ') 
w w a a 



= (-^ A)(k +k') + 



aR 



+ ( 






iR 



•Lpr 



1/2 

] B)(p+p')^^^(u +u')^^^(k+k’)^''^ 



(27) 



According to the Binomial Theorem, 



(E + e')" 

R (R - R ) 
w w a 




R -R 



w a 



(E+e’)^ 

R (R -R ) 
w w a 



(1 + 



). 



R -R 
w a 



Note that 



(28) 



(E = e')^ = + 2Ee' 



(29) 



Using Equations 28 and 29 for the left side of 27, and expanding 
the right side of 27, and neglecting small quantities fields: 



R (R -R ) 
w w a 



E r 



+ 



R (R -R )^ 
w w a 



+ 



2Ee* 



= ( 



R (R ~R ) 
w w a 



aR 



A) (k +k') + 



(30) 



But 



+ ( 



TT 2.D 



1/2 



aR 



C 1/2 

[^] B) (Puk + Puk ' 



+ pu'k + 



u k) 



(puk + puk' + pu'k + p'uk)' 






(31) 
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Thus Equation 30 becomes 



—2 E^r ' — 

, a 2Ee' 

_J_ — 



R (R -R ) R (R -R R (R -R ) 

wwa wwa wwa 



( -^ A)(k + k') + 
aR 

o 



,^1/2 Cl/2 , , , 

B)(puk)l/^lTi [f + 

o 



The time average form of Equation 32 is 



( 32 ) 




R (R -R ) 
wwa 




A)k + ( 



7t£D 

aR 



1/2 



o 



C 

_e. 

Pr 



1/2 

] B)(puk) 



1/2 



(33) 



Subtraction of Equation 33 from 32 gives 



R (R -R ) R (R -R ) o 



(34) 



wwa wwa 



+ ( 



.,1/2 C 1/2 ,—rsin , , 

,_ 2.1 .Cr.u k }. 4- li- 4- il- 



uR 



'Pr' «■ 2 



<T- + — + “^ ) • 

k u p 



The following relationships are known: 



= R at’ 



k’ = nt’ 



’ — t ’ 

P " ” ^ ~ 

a 



(35) 



Thus Equation 34 becomes 

2i 



E^ R a 
o 



R (R -R )^ 
wwa 



t’ + 



R (R -R ) 
wwa 



f / htt £ . , 

= ( — A)t’ + 
aR 

o 



( 36 ) 



+ c 



^^1/2 C 1/2 
TT i? D r_Ri D\ (t uk) 



aR Pr 
o 



B) 



(n 



t ’ t 



u 



T ) 

a u 



Since n ^ ^ << ^ 

k T u 



, Equation 36 becomes 
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r 2E 


,, . ' BP k 1 


R (R -R ) - 
i- w w a 


L 2aR u 

o 



[■ 



■ 

a R 






R (R 
w w 



1 f . 

-Rl^J 



(37) 



Equation 37 may be written 

|-R (R -R C 

w w a p 

L 4 E a R u 

o 



u* + 



R (R -R ) n7r£A 
w w a 

2 E a R 



- E R a 

o 



2 (R -R ) J 
^ w a 



t' 



(38) 



The coefficient of u* represents the sensitivity to velocity fluctuations, 
and the coefficient of t* represents the sensitivity to temperature 
fluctuations. These sensitivities could in principle be calculated, 
since all factors are either measurable quantities, physical properties, 
or known constants. In practice, however, this calculation yields values 
of the sensitivities that do not agree with those determined by actual 
calibration. This discrepancy is due to two of the supposedly known 
quantities in Equation 38, namely O' and B. King's value for B, since 
it is based on his potential flow assumption, is seriously in error. The 
value of a found in tables of physical properties may also be in error 
because the value of O' applicable to the wire in use depends very strongly 
on the history of that wire. In addition, contamination of the wire 
changes the value of D, the wire diameter, and the wire length t is some 
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effective length rather than the actual length. Thus it is impossible 
to accurately determine the sensitivities by direct calculation. In- 
stead the sensitivities are determined by calibration. 

Although the expressions for the sensitivities in Equation 38 do 
not provide quantitative information, they do indicate the direction of 
variation in the sensitivities as various parameters are changed. Thus, 
the velocity sensitivity should decrease v/ith increasing velocity, but 
velocity should not affect the temperature sensitivity to first order. 

The velocity sensitivity should increase with increasing wire temperature 
(resistance) . 

In order to determine the effect of parameter variation on the 

temperature sensitivity, further consideration of the two terms that 

comprise it is necessary. Equation 8 shows that, all other parameters 

remaining constant, an increase in T (i.e. a positive t') requires a 

a 

2 

decrease in E (i.e. a negative e') in order to maintain the equality. 
Thus, the overall sign of the temperature sensitivity must be negative. 
Since all factors are inherently positive, this implies that 



E R a R (R -R ) n Tri A 

o ^ w w a 



(39) 



2(R -R ) 
w a 



2 E a R 



Thus an increase in 



ER a 
o 



2(R -R ) 
w o 



causes an increase in 



the magnitude of the temperature sensitivity, and an increase in 
R (R -R )nirllA 

causes a decrease. 

E^ 



2E a R 



Equation 33 shows that the quantity R (R -R ) 

w w a 



is independent of T - T . That is if the parameter T - T is varied, 

w a ^ w a ’ 



/ 
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-2 - -2 

_E remains constant. Further if R is constant, E is 



R (R - R ) 
w w a 



proportional to R^. Now 



E R a 

o 



and 



w 



2 (R - R ) E 
w a 






R (R -R ) 2 

w w a 



(40) 



R (R -R ) me Ik - R (R -R ) n«M 
wwa wwa L 



2 E a R 



R cv 
o 



(41) 



Since R 



w 



and E^ 



are independent of T -T , it 

w a 



R (R -R ) 
wwa 



follows from Equation 40 that variation of T -T has no effect on 

w a 



E R a 

o 

7 /p p \ to first order. Similarly it follows from 

^ w a^ R (R -R ) n^t/A 

wwa ^ 



Equation 41 that 



2 Ea R 



increases linearly with 



increasing T -T . As previously shown, an increase in this term causes 
w a 

a decrease in temperature sensitivity. Thus an increase in T -T causes 

w a 

a decrease in temperature sensitivity. Thus it appears that velocity 
fluctuations could be measured directly by operating at very high 
temperature differences and that temperature fluctuations could be 
measured directly by operating at very low temperature differences. 
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If consideration is given to changing the dimensions and material 
of the wire, further effects on the sensitivities occur. Velocity 
sensitivity increases with wire length while the temperature sensitivity 
decreases. An increase in wire diameter causes an increase in velocity 
sensitivity and no change in temperature sensitivity to first order. An 
increase in temperature coefficient causes a decrease in velocity 
sensitivity and an increase in temperature sensitivity. 

Once the desired information on trends in sensitivity is recognized, 
it is convenient to resort to a more compact notation: 



S 

u 



4 Ea R P 



p 

1/2 - 1/2 
u 



(42) 



and 



S 

t 



ER O' 
o 



2(R -R ) 
w 



R ( R -R )n Mn 
w w a ^ 

2 Eq^ R 

o 



thus Equation 38 becomes 

e- = u- - f. 
Squaring and time averaging yields 



(43) 



(44) 



2 2 2 2 

e’ = S - 2S u’t’ + S t' . 

u u T T 

The sensitivities S and must be determined 

u T 

This is discussed in the section on procedure. The 



(45) 

from calibration. 



quantity, e’ , 
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is experimentally measured. From known values of e * ’ ^ 

of simultaneous equations of the form of Equation 39 can be formed. 

2 2 

This set can then be solved for u* , u*t*, and t' . This solution is 
discussed in the RESULTS section. 



1 




I 



I 
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APPARATUS 



Flow System 

The flow system is shown schematically in Figure 1. A Craftsman 
Model 315« 16970, 3HP industrial vacuum cleaner was operated as a blower 
to supply air to the system. A three inch galvanized duct carried the 
air to the highest point of the flow system where the duct tapered to 
meet the 1.0625 inch inside diameter Type M copper tubing that comprised 
the remainder of the loop. The three- inch duct was fitted with a bleed 
line and damper. The flow rate in the loop could be varied by using the 
damper to alter the division of air between loop and bleed line. Above 
the bleed line in the duct was an orifice meter for measuring air flow 
rate in the loop. (The calibration curve for the orifice meter is 
given in Appendix B.) In the copper tube, the air passed downward 
through an 83 inch unheated velocity developing section. This was 
followed by a 58 inch heating section which was wound with nichrome 
heating wire and covered with 1.125 inches of pre-formed pipe insulation 
Probes were inserted upward into the copper tube as shown. All measure- 
ments were made in the end plane of the heating section. 

Ir on-constantan bulk thermocouples were installed at the beginning 
of the velocity developing section and at the end of the heating section 
and wall thermocouples were installed at the end of the heating section. 



I 

I 
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Figure 1. Flow Diagram of Apparatus 
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Probes and Traversing Mechanisms 

The rapid response thermocouple used (a Heat Technology Laboratory 
TCFW 202-ChA-5) is shown in Figure 2. The response time of this thermo- 
couple has been estimated by Rodr iguez-Ramirez to be 5 milli-seconds . 

The hot-wire probe used (a Thermo Systems 1210-T1.5) is shown in 
Figure 3. 

Positioning of the rapid response thermocouple was by means of the 
traversing mechanism shown in Figure 4. When in place as shown in Figure 
1 the mechanism was, of course, inverted. The ball joint permitted 
radial movement of the probe across the tube. Radial position was 
measured with the dial guage shown. 

Figure 5 shows the traversing mechanism used to position the hot- 
wire in the loop. This traversing mechanism replaced the thermocouple 
traversing mechanism at the bottom of the loop when hot-wire measure- 
ments were being made. The upper block in Figure 5 could be moved 
laterally so as to place the hot-wire probe at the desired radial position, 
while the lower plate shown in Figure 5 was bolted to the loop. A pitot 
tube and a local mean thermocouple were inserted through the auxiliary 
probe entrances. 
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Figure 2. Fast Response Thermocouple 
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Figure 3. The Hot-Wire Probe 
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Figure 4. Thermocouple Traversing Mechanism 
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Figure 5. Hot-Wire Traversing Mechanism 
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Instrumentation 

Instrumentation for hot-wire measurements consisted of a Thermo 
Systems Type 1010 anemometer, and a Ballantine Model 320A True Root- 
Mean-Square Electronic Voltmeter, 

For rapid- response thermocouple measurements, instrumentation consist- 
ed of a Tektronix 122 Amplifier (See Appendix C for frequency response.) 
and the Ballantine Voltmeter. 

Bulk and wall thermocouple measurements were made with a Dana Model 
2850 D.C. amplifier and a Beckman Model 4011R digital voltmeter. 
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PROCEDURE 



Hot-wire in Isothermal Flow 

The first step was calibration. The constant C defined by Equation 
17 and appearing in Equation 18 can be determined from mean voltage and 



were plotted as a function of u^, the result should be a straight line 
with slope C . 

The Calibration of the hot-wire probe was carried out with the probe 
in place in the loop and located at the tube centerline. Variation of 
u was achieved by varying the flow rate in the loop using the bleed line 
damper. The quantity E was measured at seven or eight different flow 
rates. The value of u at the centerline corresponding to each flow rate 
was measured with a pitot tube. The data were plotted as indicated above, 
a straight line was fitted to the points, and the slope was measured. 
Calibration curves are shown in Appendix E. 

The second step was the making of fluctuation measurements. Equation 
18 gives the relationship between the rms velocity fluctuation and the 
rms voltage fluctuation. 

The flow system was operated at Reynolds Numbers of 30,700, 41,500, 



and 49,800. For each Reynolds Number, measurements were made of\e' 



velocity measurements. It is seen from Equation 14 that if E 




2 



at eight points between the center and the wall of the tube corresponding 
to r/R of 0, 0.125, 0.250, 0.375, 0.500, 0.625, 0.750 and 0.875. A 
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pitot tube was used to measure u at the same points. 

After all fluctuation measurements had been made, the calibration 
procedure was repeated to make certain that the constant C had not 
changed during the run. 

Rapid Response Thermocouple in Non- Isothermal Flow 

The flow was adjusted so that the Reynolds Numbers in this part were 
the same as those in the previous section. The applied heating was 430 
watts in all cases. 



Measurement was made of from the rapid response thermocouple 

as a function of radial position for each Reynolds Number. No calibration 

was done, since the constant, a, in Equation 24 was assumed to be that 

13 

found in standard thermocouple tables . This assumption is of question- 
able validity. 

A conventional thermocouple was used to measure the centerline 
temperature, and the wall thermocouples were used to measure the wall 
temperature for each Reynolds Number. 



The procedure of calibration and measurement outlined below repre- 
sents one of the contributions of this work. It is anticipated that 
this procedure will be used for the glycol and mercury measurements 
previously mentioned. 

As shown in the THEORETICAL BACKGROUND section, the sensitivities 

S and depend on two parameters: u and T -T . Thus, in calibration 

u T w a ’ 

it is necessary to determine the functional dependence on each parameter. 
Calibration for temperature sensitivity was also carried out with 
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the probe in place and at the tube centerline. Velocity variation was 

again achieved by use of the bleeder line damper, and u was measured with 

a pitot tube. In this way five different values of u were used, and 

for each value of u thirteen values of T -T were used. 

w a 

These data were plotted as E versus u for constant values of T -T 

w a 

and as E versus T -T for constant values of u. These curves are shown 
w a 

in Appendix G. The sensitivities and are the derivatives of these 

curves. These derivatives were obtained from analytical differentiation 

of parabolas fitted to the data by the method of least squares. The 

selection of parabolas is arbitrary, but the relatively slow and smooth 

variation of the data makes a parabola a satisfactory approximation. 

Note that T -T must be computed from measurements made of R -R . 

w a ^ w a 

The resistance- temperature characteristic of the hot-wire probe was 

experimentally determined to be as shown in Appendix D. 

Once S and were known as functions of u and T - T , fluctuation 
u T w a 

measurements were made. The rms voltage fluctuation was measured at 
r/R = 0 and r/R = .5 for each flow rate and temperature difference con- 
sidered in calibration. 

When all fluctuation measurements had been completed, the calibration 
was repeated to determine the drift in and . 
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RESULTS 



Hot-Wire in Isothermal Flow 

Figure 6 shows the axial turbulent intensity at the centerline of 
the pipe compared with the work of Sandborn. 

Figure 7 shows axial turbulent intensity as a function of radial 
position for three different Reynolds Numbers compared with the work of 

T ^ 

Laufer . 

The agreement of present results with those of Sandborn^^ and 
Laufer^^ is rather good. The data used by Sandborn to obtain his 
correlation show a scatter of t 10 percent. The results of the present 
work are well within this range. 



Rapid Response Thermocouple in 
Non- Isothermal F low 

Figure 8 shows the temperature fluctuation intensity as a function 

14 

of radial position. Rodr iguez-Ramirez , and also Tanimoto and Hanratty 
have published results of this nature. The present results are for a 
range of Reynolds Numbers higher than the range considered by Rodriguez- 
Ramirez; however, a preliminary run was made in which the flow rate and 
heat input were the same as those in one of Rodr iguez-Ramirez * 
runs. The results of this run deviated from those of Rodr iguez-Ramirez 
by less than 5 percent. The temperature gradients used by Tanimoto and 
Hanratty were so much larger than those used in the present work that 



comparison is not possible. 
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F i gure 7 . 



Velocity Turbulence Intensity in Air I 



29 




Figure 8. Temperature Turbulence Intensity in Air 
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Hot-Wire in Non- Isothermal Flow 

Figures 9 and 10 show the dependence of the sensitivities S and 

on the parameters u and T -T . The velocity sensitivity, S , increases 

w a u 

with increasing temperature difference and decreases with increasing 
velocity as predicted from Equation 38. The temperature sensitivity 
decreases with increasing temperature difference as predicted from Equation 
38. Note that there is a secondary dependence on velocity not predicted 
by Equation 38. This dependence is due to the velocity dependence of the 
effective wire length. 

The investigation of various calculating schemes including the 

selection of the best one for analyzing the non- isothermal data was an 

important part of the present work. As pointed out in the section on 

theoretical development, Equation 45 gives the relationship between the 

2 2 

unknown quantities u’ , u't’ and t' and the known quantities S , S^, 

2 

and e' . Since there are three unknowns, at least three equations of the 

form of Equation 39 are needed. These can be obtained by operating the 

hot-wire at least three values of T (and hence of T -T ) . For each 

w w a 

value of T -T , the values of S and S appropriate to the mean velocity, 
w a u T 

2 

u, of the flow are read from Figures 9 and 10, and e’ is measured. Thus, 

for each value of T -T an equation is formed. 

w a 

Once the set of equations is obtained, a number of methods of solution 

are available. The most obvious one is the simultaneous solution of 

16 

three such equations by a method such as Gaussian reduction. However, 
the numerical calculations involve the subtraction of large numbers whose 
difference is small, which has the effect of greatly magnifying the errors 

in the data. Since the errors in and were already large, the result 

^ 2 

, u't’, and t ' of the 



was errors in the calculated values of u' 
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Figure 9 



Hot-Wire Temperature Sensitivity 
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F igure 10 



Hot-Wire Velocity Sensitivity 
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order of 500 percent . 
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The Kovasznay Fluctuation Diagram is a graphical method of solving 

2 

the set of equations in which equation 39 is divided by to give: 






(46) 



p.2 S 

e u 

Then the quantity r is plotted versus ■— . Note that analytically 

St t 

the relationship between the two is described by a parabola. By the 
method of least squares a parabola is fitted to the data. The computer 
program used to do this fitting is listed in Appendix F. The intercept 



,2 



of this parabola with the 



— axis is t^^ , the slope at -^=0is 

Z 



:>T 2 ^ 

“ 2u*t^ , and the curvature of the parabola is u^ . This method was 

S 

used on the present data, but the need to extrapolate the data from ■— of 

S S ^T 

a u 

some non-zero value to - =0 and then determine the slope at ■— = 0 

resulted in a lar<^e uncertainty in u’t’ . This same problem was encountered 
9 

by Arya and Plate , who suggested that if an independent measurement of 

2 

t* could be obtained to locate the intercept of the curve accurately, 



,.2 



then reliable values of the remaining two unknowns u’t^ and u*^ could 
be determined by the Kovasznay method. This procedure was followed with 



.2 



the present data. The value of t used was the value found using tl^e 
rapid response thermoc oup le . The results were still rather uncertain, 
but the uncertainty was less than for any of the other methods. (The 
Kovasznay Fluctuation Diagrams for representative cases are shown in 
Appendix H.) 



2 ~ 

The numerical values of f u’ t presented in Figures 11 and 12 
deviate seriously from the isothermal results wIlii wnicli it is thought 
they should agree. (See, for instance Corrsin and Uberoi ) The trends 
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of the results with Reynolds Number and radial position are in accordance 
with expectations. 

The calculated values of u ' t ' /u^ (T^-T^) are presented in Figure 13. 

Turbulence at the centerline in pipe flow for the Reynolds numbers used 

in the present work should be isotropic. Thus u't' and u ' t ' /u (T -T ) 

L oh 

should both be zero at the centerline. The centerline results in 
Figure 13 do not exhibit this behavior. However the error in these re- 
sults is roughly estimated in Appendix I as 1007c. 

Away from the centerline the turbulence should no longer be iso- 
tropic, and u't' should no longer be zero. With heating from the wall 

as in the present case, u't', and hence u't'/u^ (T -T ) should be negative. 

cot 

This trend is followed by the present results in that u't'/u (T -T ) at 

tot 

r/R = .5 has a larger magnitude than at r/R = 0 and in that it has a 

negative sign. The quantity u't’/u, (T -T^) cannot continue to increase 

tot 

indefinitely as Reynolds Number increases, for this would mean no dis- 
sipation was taking place. The increasing trend in the results shown in 
Figure 13 for r/R = .5 is therefore not anticipated, but it is possible 

so long as the trend reverses itself at higher Reynolds Numbers. Arya 
9 

and Plate found that in the boundary layer growing on the wall of a wind 
tunnel u' t ' /u^(T^-T^)remained nearly constant with increasing Reynolds 
Number. Again with an estimated error of 100 percent little significance 
can be attached to present results. 
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Figure 13 . 



Temperature-Velocity Turbulence Intensity in Air 



f 

i 



f 



38 



DISCUSSION 

The numerical results forV u’^/u^ and u't*/u, 

t L 

are in general, too large. In view of the way in which these intensities 

are obtained from the Kovasznay Fluctuation Diagram, this suggests that 

S is too small. The temperature sensitivity affects both coordinates 
u T 

in the same way, and therefore, errors in S would be self-compensating. 

S ^ 

However, errors in S affect only the — coordinate. Too small a value 

T 

of S would cause too large a slope and curvature in the fitted parabola 

2 

and too large a value of u' and u*t'. 

Because of the long time required for the system to come to thermal 

equilibrium after each change in flow rate, the data used to compute 

for a given value of T -T were taken over a period of approximately 

w a 

eight hours. It is reasonable to expect an increase in the thickness of 

the layer of contamination on the wire during this period since the wire 

was continuously exposed to an inadequately filtered flow. 

3 

Kronauer has shown that dirt accumulations on fine wires at low 
flow rates increase the heat transfer rate because the heat transfer area 
increases more rapidly than the heat transfer resistance as dire 
accumulates. Thus, a dirty wire would require a higher value of E than 
a clean one to maintain it at a given temperature. 
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The wire used for the present measurements was examined under a 
microscope, and a significant amount of dirt was observed on it. 

In the present calibration, E was measured at the highest calibration 
velocity first, and the velocity was decreased step by step until the 
lowest value was reached. As anticipated, E decreased with decreasing 
velocity, but the increasing layer of contamination apparently worked 
counter to this trend. Thus, in the calibration procedure used the rate 
of decrease of E was slower than if contamination had not been a problem. 
The result was a value of S that was too small. 

Larger values of S would result in larger values of at the same 

2 2 ^ 

value of e * - t' on the Kovasznay Fluctuation Diagram and thus, a 




smaller slope and curvature of the fitted parabola. The calculated values 



It was not possible to determine the errors in and u’t' 

caused by the drift in E noted above because it was not possible to 
measure what E would have been without contamination. Thus it was 
necessary to use a different method to estimate these errors. The method 
used is outlined in Appendix 1, and the error estimates are 36 percent 
in V 100 percent in u't* . As noted in Appendix 1, these figures 

do not represent true confidence limits on the results. They are only 
rough estimates of the errors present. 

The calibration curves for before and after as shown in Figures 24 
through 26 show such a large drift that all data taken in between should 




and u't' would then be smaller and in better agreement with the 



isothermal results. 




be discarded. This was not done, since a number of attempts to obtain 
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better results failed. Note then, that all fluctuation intensities ob- 
tained from the present non- isothermal hot-wire measurements are of very 
doubtfu 1 reliability. 

Due to inexperience in this area, a number of deviations from good 
experimental practice were present. First, although the wire was 
examined under a microscope after all measurements had been made, no 
examination was made before or during experimental runs. Thus the 
hypothesis of an increasing dirt layer cannot be verified. Second, no 
attempt was made to clean the wire during the course of experimental 
runs. Because of these errors it is not possible to analyze in detail 
or to correct for the drift that occurred. 

In future work, provision should be made for careful filtration of 
the fluid and a means of rapid calibration to minimize drift. Two cali- 
brations, one with increasing velocity and one with decreasing velocity, 
should be made to check for contaminat ion . If it is still not possible 
to reduce drift to an acceptable level, provision should be made for the 
complete withdrawal of the probe from the flow. 

A problem encountered with the present apparatus was difficulty in 
maintaining flow rate and heat input at desired levels. Long term 
temperature variations of two or three degrees persisted even after the 
system had come to ’’equilibrium". This was apparently caused by cycling 
of the room temperature. 

The nature of the problems encountered in the present work serve as 
a guide in planning additional steps in the overall research program. 



41 



SUMMARY AND CONCLUSIONS 

This work has developed procedures for the calibration of a hot- 
wire in non-isothermal flow. Further, the various calculation schemes 

available for analyzing non-isothermal flow fluctuation data have been 

9 

evaluated and the Arya and Plate modification of the Kovasznay Fluctuation 
Diagram^^ method has been found to be the most satisfactory. 

Apparatus limitations and deviations from good experimental practice 
precluded the production of reliable numerical results. 

The findings of this investigation are an important preliminary step 
in a comprehensive research program of non-isothermal flow measurements 
with hot-wire and hot- film anemometers. 
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APPENDIX A 



APPENDIX A 



TABLES OF EXPERIMENTAL DATA 



Table I 

Orifice Meter Calibration 

AP Orifice AP Pitot Tube 

Meter (in)* At (in)* 



4 . S53 


2.173 


3.151 


1.595 


2.172 


1 . 123 


1.114 


0.571 



* Manometer Fluid SPC. 0.797 



I 
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Table 2 

Tektronix Amplifier Frequency Response 

Frequency Voltage in Voltage Out 

(CPS) (Millivolts) (Volts) 



5 


.200 


.200 


7 


.200 


.214 


9 


.200 


.222 


20 


.200 


.238 


40 


.200 


.244 


60 


.200 


.245 


80 


.200 


.245 


100 


.200 


.244 


200 


.200 


.239 


300 


.200 


2.31 


400 


.200 


.226 


500 


.200 


.217 


600 


.200 


.210 


700 


.200 


.200 


800 


.200 


. 193 


900 


.200 


. 185 


1000 


.200 


. 177 


1100 


.200 


. 170 
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Table 3 



Calibration for Isothermal Re=49,800 Run 



P ( in of 

.tot Tube SPG. 797 
F luid ) 


R (OHM) 
a 


R (OHM) 
w 


E (VOLTS ) 


/ e ' ^(millivolts 


3.890 


7.57 


10.60 


3.07 


13.0 


3.394 


7.54 


10.56 


3.04 


13.0 


2.948 


7.52 


10.53 


2.98 


12.9 


2.407 


7.50 


10.50 


2.93 


12,8 


1.935 


7.48 


10.47 


2.88 


12.7 


1.369 


7.46 


10.44 


2,81 


12.5 


0.844 


7.43 


10.40 


2.71 


12.5 


0.330 


7.38 


10.33 


2.54 


12.6 


fore * 










4.004 


7.57 


10.60 


3.08 


— 


0.737 


7.42 


10.39 


2,70 


~ ~ . 



After 



t 
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Table 4 





Data for 


Isothermal Re=49,800 Run 




Radial 


ap 


(in of 


o 


Position (in) 


Pitot 


Tube SPG. 797 E(Volts) 


\ e' (Millivolts) 



Fluid) 



0.0000 


3.984 


3.07 


12.9 


0.0625 


3.964 


3.07 


12.9 


0. 1250 


3.841 


3.07 


16.5 


0.1875 


3.691 


3.06 


19.5 


0.2500 


3.433 


3.04 


23.6 


0.3125 


3.184 


3.02 


27.5 


0.3750 


2.870 


2.97 


31.7 


0.4375 


2.411 


2.94 


36.4 



Table 5 



Calibration for Isothermal Re=30,700 and 41,500 Runs 



AP Pitot Tube 

(in. of SPG. 797 Fluid) 


R^(Ohms) 


R (Ohms) 
w 


i: (Volts) 


2 

f e' (Millivolt 


3.533 


7.55 


10.57 


3.05 


12.8 


2.999 


7.53 


10.54 


2.98 


12.6 


2.500 


7.51 


10.51 


2.94 


12.5 


2.000 


7.48 


10.47 


2.88 


12.4 


1.496 


7.46 


10.44 


2.82 


12.4 


1.020 


7.44 


10.42 


2.75 


12.3 


0.500 
Before ^ 


7.40 


10.36 


2.61 


12.4 


2.780 


7.52 


10.53 


2.96 


— 


1.020 
After . 


7.44 


10.42 


2.75 
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Table 6 



Data for Isothermal Re=30,700 Runs 



Radial .. . 

Position 


AP Pitot Tube 
(in. SPG. 797 Fluid) 


E (Volts) 


^ e ' ^ (Mi 1 1 ivolt s ) 


0.0000 


1.480 


2.81 


12.2 


0.0625 


1.460 


2.81 


12.9 


0.1250 


1.400 


2.80 


15.1 


0.1875 


1.330 


2.79 


18.2 


0.2500 


1.240 


2.78 


21.9 


0.3125 


1.130 


2.76 


24.8 


0.3750 


1.020 


2 .74 


28.6 


0.4375 


0.850 


2.71 


32.5 



Table 7 

Data for Isothermal Re=4 1,500 Run 



Radial . 4P Pitot Tube 

Position (in. SPG797 Fluid) E (Volts) . e’ (Millivolts) 



0.0000 


2.780 


2.95 


12.4 


0.0625 


2.730 


2.95 


13.4 


0.1250 


2.630 


2.95 


15.6 


0.1875 


2.500 


2.94 


19.0 


0.2500 


2.330 


2.92 


22.6 


0.3125 


2.150 


2.90 


26.2 


0.3750 


1.960 


2.89 


29.4 


0.4375 


1 .650 


2.85 


34.2 
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Table 8 

Thermocouple Data for Re=30,700 Run 



Radial 

Position 

fin) 



0.00 



0.05 



0.10 0.15 



0.20 0.25 



0.30 0.35 



0.40 0.45 NOISE 



/ e'2 

(Mi 1 li- 

volts) 20.0 21.0 23.5 28.5 30.5 34.5 38.5 42.5 45.0 47.5 8.0 



Table 9 

Thermocouple Data for Re=4 1,500 Run 

Radial 

Position 

(in) 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 NOISE 



(Milli 

volts) 15.0 15.7 18.3 21.5 25.0 28.5 32.0 35.5 38.0 39.5 6.0 



Table 10 

Thermocouple Data for Re=49,800 Run 



Radial 

Position 

(in) 0.000.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 NOISE 



V e ' 2 

(Milli- 
volts) 12.0 12.8 14.8 17.5 20.5 24.0 27.0 30.0 32.0 33.2 3.0 
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Table 11 

Non- Isothermal Hot-Wire"Before"Calibrat ion Data 
AP (in)"^ 

R -R (Ohm) 2.20 1.620 1.330 0.932 0.480 



a i 




1.00 


1 


.865 


1.820 


1 


.795 


1. 


,740 


1 


. 665 


1.20 


1 


.990 


1.940 


1 


.915 


1, 


,850 


1 


.775 


1.40 


2 


125 


2.070 


2 


.045 


1. 


,955 


1 


.875 


1.60 


2 


.230 


2.175 


2 


. 145 


2. 


,080 


1 


.960 


1.80 


2 


.325 


2.265 


2 


.235 


2. 


, 180 


2 


.075 


2.00 


2 


.410 


2.355 


2 


.320 


2. 


.265 


2 


. 150 


2.20 


2 


. 500 


2.435 


2 


.400 


2. 


,340 


2 


.225 


2.40 


2 


. 575 


2.515 


2 


.475 


2, 


.410 


2 


.295 


2.60 


2 


.650 


2.590 


2 


.550 


2, 


.490 


2 


.360 


2.80 


2 


.725 


2.660 


2 


.620 


2 


.555 


2 


.425 


3.00 


2 


.790 


2.725 


2 


.680 


2 , 


.620 


2 


.485 


3.20 


2 


.855 


2.785 


2 


.745 


2 


.686 


2 


. 585 


3.40 


2 


.920 


2.84 5 


2 


.805 


2 


. 746 


2 


.640 



Orifice Meter Manometer Fluid SPG 2.95 
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Table 12 

Cold Resistances in" Bef ore"Calibration 



‘ (in.)* 


R ..(OHM) 
cold 


2.120 


7 .42 


1.620 


7.46 


1.330 


7.49 


0,932 


7.57 


0.480 


7.80 
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Table 13 



Temperature Difference Between Wall and Centerline 



AP (in)* 



T - T ”F 

r/R=l r/R=0 



2.120 

1.620 

1.330 

.932 

.480 



46 

51 

54 

61 

74 



*Orifice Meter Manometer Fluid SPG 2.95 
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Table 14 

Non-Isothexmal Hot-Wire RMS Voltage 
Fluctuation Data for r/R = 0 



AP (in)*-> 



R -F (Ohm) 
w a 

1 


2. 120 




1.620 


1.330 


0.932 


0.480 


1.00 


24.8 


Milli- 

volts 


27.0 


28.8 


31.6 


38.3 


1.20 


24.8 




26.9 


28.4 


31.0 


37.2 


1.40 


24.7 




26.3 


28.2 


30.3 


36.2 


1 .60 


24.5 




26.0 


27.9 


29.8 


35.3 


1.80 


24.3 




25.7 


27.4 


29.2 


34.3 


2.00 


24.3 




25.5 


26.9 


28.9 


33.6 


2.20 


24.2 




25.2 


26.7 


28.5 


32.9 


2.40 


24. 1 




25.0 


26.4 


28.1 


32.0 


2.60 


24.1 




24 .9 


26.2 


27.7 


31.4 


2.80 


24.0 




24.7 


25.9 


27.5 


31.0 


3.00 


24.0 




24.5 


25.8 


27.2 


30.6 


3.20 


23.9 




24.5 


25.6 


27.0 


30.0 


3.40 


23.9 




24 .4 


25.5 


26.7 


29.8 
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Table 15 



Non- I s otherma 1 Hot-Wire RMS Voltage 
Fluctuation Data for r/R=.5 

AP(in)* > 



R -R (Ohm) 
w a 


2.120 




1.620 


1.330 


0.932 


0.480 


1.00 


42.9 


Milli- 

volts 


45.0 


46.8 


50.5 


59.7 


1.20 


42.9 




44.8 


46.8 


50.0 


58.3 


1 .40 


42.6 




44.6 


46.2 


49.2 


57.1 


1.60 


42.6 




44.0 


45.8 


48.7 


56.2 


1.80 


42.7 




44.0 


45.6 


48.0 


55.0 


2.00 


42.7 




43.8 


45.1 


47.5 


54.2 


2.20 


42.3 




43.5 


44.9 


47.0 


52.8 


2.40 


42.3 




43.3 


44.8 


46.2 


52.2 


2.60 


42.2 




43.2 


44.5 


46.0 


51.7 


2.80 


42.8 




43.2 


44.5 


45.8 


50.8 


3.00 


42.6 




43.0 


44.3 


45.2 


50.3 


3.20 


42.6 




43.0 


44.3 


45.1 


49.5 


3.40 


42.6 




43.0 


44.0 


44.8 


49.2 



* Orifice Meter Manometer Fluid SPG 2.95 
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Table 16 

Non- Isothermal Hot -Wirc’^'Af ter"Ca 1 Ibrat Ion Data 
AP (in)* ^ 



R -R (Ohm) 
w a 


2.120 


1.620 


1.330 


0.932 


0.480 


1.00 


2.010 (Volts) 


1.945 


1.920 


1.860 


1.775 


1.20 


2. 140 


2.100 


2.070 


1.985 


1.890 


1.40 


2.235 


2.215 


2.180 


2.120 


2.000 


1.60 


2.365 


2.320 


2.285 


2.220 


2.125 


1.80 


2.470 


2.420 


2.385 


2.320 


2.215 


2.00 


2.565 


2.515 


2.480 


2.405 


2.295 


2.20 


2.655 


2.600 


2.565 


2.495 


2.375 


2.40 


2.740 


2.685 


2.645 


2.570 


2.450 


2.60 


2.820 


2.765 


2.720 


2.645 


2.525 


2.80 


2.900 


2.840 


2.795 


2.720 


2.595 


3.00 


2.970 


2.915 


2.860 


2.785 


2.655 


3.20 


3.070 


2.980 


2.930 


2.845 


2.715 


3.40 


3.135 


3.075 


2.995 


2.910 


2.775 
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Table 17 

Cold Resistances in”Af ter’*Calibrat ion 



' (in.)* 


R (OHM) 

cold 


2.120 


7.21 


1.620 


7.24 


1.330 


7.26 


0.932 


7.29 


0.480 


7.43 



,1 






1 

I 

i 



I 






i 



I 
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APPENDIX B 

ORIFICE METER CALIBRATION 

The orifice meter used was calibrated in place against a pitot tube 

located at the centerline of the test section. von Karman ’ s universal 

1 8 

velocity distribution was used to calculate flow rates from the pitot 

tube velocity measurements. The calibration curve is shown in Figure 14. 

It is found as anticipated that a plot of flow rate versus the square 

19 

root of the pressure drop is a straight line. 






I 



FLOW RATE ( LB / HR ) 
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Figure 14. 



Orifice Meter Calibration Curve 
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APPENDIX C 



AMPLIFIER FREQUENCY RESPONSE 

The frequency response of the Tektronix amplifier is shown on 
the following page. 



I 



1 
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Figure 13. Amplifier Frequency Response 
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APPENDIX D 

RESISTANCE -TEMPERATURE CURVE FOR THE HOT-WIRE 

Figure 16 shows the relationship between resistance and temperature 
for the hot wire used in the non- isotherma 1 flow measurements. The 

3 

dashed curve is the approximate relationship suggested by Kronauer. 



7 
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APPENDIX E 

CALIBRATION OF HOT-WIRE IN ISOTHERMAL FLOW 

The isothermal flow data were taken on two different days. The set 

of data for a Reynolds Number of 49,800 was taken on the first day. The 

calibration curve for these data is shown in Figure 17. The constant C 

-2 

in Equation 18, which is the slope of the line in Figure 17 is 1.67x10 

The data at Reynolds Numbers of 30,700 and 41,500 were taken on the 

second day. The calibration curve for these data is shown in Figure 18. 

-2 

The constant C for these data is 1.69x10 

-2 

The value of C calculated from Equation 17 is 0.84 x 10 . The 

discrepancy between this value and that obtained from calibration is due 

to wire historv effects in and R and the error Inherent in King's 

o 

potential flow value for B. 



( VOLT / OHM ) 
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Figure 17. Isothermal Calibration of Hot-Wire I 
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Figure 18. Isothermal Calibration of Hot-Wire II 
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APPENDIX F 

LISTING OF CURVE FITTING PROGRAM 

The computer subroutine used to perform least square curve fitting 
calculations is listed on the following pages. 
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APPENDIX G 

CALIBRATION OF HOT-WIRE IN NON- ISOTHERMAL FLOW 

Figures 19 through 23 represent the calibration data used to compute 
the results presented for non- isothermal flow. The curves drawn through 
the data represent the least squares fit of parabolas to the data. These 
parabolas were differentiated analytically to obtain the values of and 
presented in Figures 9 and 10. 

Figures 24 through 26 show calibration curves for representative 

values of u and T - T before and after fluctuation measurements were 
w a 

made. Table 18 gives a complete tabulation of and as functions of 

u and T - T before and after the fluctuation measurements, 
w a 
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Figure 19. Non- Isothermal Calibration of Hot-Wire I 



E ( VOLTS ) 
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u ( FT / SEC ) 



Figure 20. Non- Isothermal Calibration of Hot-Wire II 
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lUJ 



F igure 




21. Xon- Isot hernial CalibraLion of Hot-Wire III 



( VOLTS ) 
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Figure 22. Non- Isothermal Calibration of Hot-Wire IV 
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Figure 23. Non-Tsothermal Calibration of Hot-Wire V 



( VOLTS ) 



76 




Figure 24. Calibration Drift I 
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Figure 25. Calibration Drift II 
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igure 26. Calibration Drift III 
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APPENDIX H 



kovasznay fluctuation diagrams 



Figures 27 through 29 show the Kovasznay Fluctuation Diagrams for 
three representative cases. 

The parabolas that fit these data are: 

F igure Parabola 



.2 



S 2 



27 



28 



29 



( - t'^) = 2.64 +16. 9(^) + 19.1(^) 

^T ^T 



. 2 , 



S 2 



( — R - t'") = 0.43 +18. 9(-^) + 16. 2(^) 

^T T 



2 S S 2 

( - t’^) = 0.21 + 42. 6(-^) + 47. 3(^) . 

^T ^T 



u' is then found by dividing the coefficient of ~ by -2, 

— 2 ^ 

and V u* is the square root of the coefficient of ( ~) . The 



results from the parabolas given above are listed in Table 19. 



Table 19 



Sample Results for v u’ and u’t* 



Reynold s 
And r/R 


Number 


. u-2 


u* t * 


27,400 


0.0 


4.36 


-8.45 


38,200 


0.0 


4.03 


-9.45 


38,200 


0. 5 


6.89 


-21.3 
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Figure 28. 



Kovasznay Fluctuation Diagram II 



i1 
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Figure 29. 



Kovasznay Fluctuation Diagram III 
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APPENDIX I 



ERROR ANALYSIS 



The results of the isothermal flow measurements with the hot-wire 
deviate from the results of Sandborn^^ and Laufer^^ by less than 10 
percent . 

The results of a preliminary non- isothermal run duplicating the 
conditions of Rodr iguez-Ramirez and using the rapid response thermocouple 
agreed with his results within 5 percent. 

Since the errors in the results of the non- isothermal flow measure- 
ments with the hot-wire are quite large, a detailed estimate of their 
magnitude is needed. 



crepancies between successive calibrations as tabulated in Table 18. 

Note that the actual errors in S and S_ come from a continuous 

u T 

drift in E during a given calibration as outlined in the DISCUSSION 
section. Since it is not possible to measure this drift, recourse was 
made to analyzing by this method. Thus the error estimates obtained 
below do not represent . confidence limits on the results but only rough 



It is assumed that the errors 




are primarily the result of errors in S and S^ . In order to obtain an 
^ u T 

estimate of the errors in S and S,^ , an analysis was made of the dis- 

u r 



approximations to the confidence limits. In Table 18 , the largest 



86 



discrepancies occur for u = 137 ft/sec. The percentage error, based on 
the before value, between the before and after values, was determined 
for each entry in the 137 ft/sec column. All the errors were summed 
and divided by the number of entries and likewise for . Thus, average 
errors were determined to be 19 percent in and 14 percent in . 

The propagation of these errors into the coordinates of the Kovasznay 
Fluctuation Diagram is as follows: 



A ( e 



.2 






,2 



AS, 



= 28 percent , 



A (iu) 



^ ^u + ^ ^T = 33 percent. 



These errors were assumed to apply to all values of u and r/R. 

S 

The '^alue of for each point of the Kovasznay Fluctuation 

^T 

Diagram for u = 66 ft /sec r/R = 0 was increased by 33 percent, and the 



2 2 

value of e * - t' was decreased by 28 percent. A parabola was fitted 



to the resulting points and u’ and u‘t' determined. The new values 
differed from the original values by 72 percent and 100 percent respec- 
tively. Since 



V u'^ = I . u' 



72 2 ^ 



i 
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the error in V ^ ^ is 36 percent. This 36 percent error was used as 



the limit of error for all values of \u* /u^ , and 100 percent was used 

t 

for all values of u't'/u _ (T ) . 

cL o L 



